Model for gravitational collapse in effective quantum gravity 
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We describe a non-perturbative approach to studying the gravitational collapse of a scalar field 
in spherical symmetry with quantum gravity corrections. Quantum effects are described by a phase 
space function that modifies the constraints of general relativity but does not affect the closure of 
their algebra. The model may be viewed as one that incorporates a class of anomaly-free quantum 
gravity effects. Numerical simulations of the resulting field equations show that the model repro- 
duces known classical results for sufficiently massive initial data but gives a dramatic change at the 
threshold of black hole formation characterized by oscillations and a mass gap. 
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General relativity (GR) is incomplete as a theory of 
gravity. This is because the gravitational collapse of 
physical matter produces a black hole with a singularity 
inside, which is a point in spacetime where the curvature 
and energy density diverge. The hope and expectation is 
that this problem will be overcome in a quantum theory 
of gravity. Evidence for this comes from both physical 
arguments and the quantization of cosmological models 
with singularities. 

Gravitational collapse in classical GR is well studied 
numerically in the spherically symmetric case with a va- 
riety of matter fields. The first such model was with 
a massless minimally coupled scalar field [j,]. This pio- 
neering work provided a mass scaling law at the onset 
of black hole formation which indicated that black holes 
form with infinitesimal mass. It also provided a finely 
tuned violation of the cosmic censorship conjecture, and 
a discrete self-similarity of the field variables that extends 
to arbitrary scales. 

These results are not expected to hold in a quan- 
tum theory of gravity because of quantum fluctuations 
of the metric in regions of sufficiently large energy den- 
sity. The general understanding derived from homoge- 
neous isotropic models is that singularity avoidance is 
associated with violation of energy conditions, which in 
turn suggests that quantum gravity provides a repulsive 
force at very short distances. 

These considerations suggest that it should be of much 
interest to find the quantum theory of GR coupled to 
the massless scalar field theory in spherical symmetry. 
However this model is a non-trivial two-dimensional field 
theory for which no quantization is known, in either the 
canonical quantization program or in the string theory 
approaches to quantum gravity. Given this state of af- 
fairs, we introduce and study in this letter a model that 
modifies this sector of GR by taking into account a pos- 
sible class of quantum effects. 

The class of modifications we consider is motivated by 
the effective constraints approach used in the loop quan- 
tum gravity (LQG) program in the context of cosmology 
There are two general types of quantum corrections 
in LQG, which come from the inverse triad and holon- 
omy operators [3|. The work in cosmology attempted 



to incorporate the inverse triad corrections, and has re- 
cently been extended to the spherically symmetric sector 
in LQG 

Li this letter we introduce a similar effective approach 
in the Arnowitt-Deser-Misner (ADM) canonical formu- 
lation. This is motivated by the polymer quantization 
method applied to the spherically symmetric sector of 
GR 0, 01 ■ In this formulation there is an analog of the 
inverse triad operator of LQG, which is the operator re- 
alization of the factor I / ^ det {[h])' in the Hamiltonian 
constraint of GR (where h is the spatial metric). 

Effective quantum gravity corrected equations based 
on this observation are investigated in @ and where 
the corrections are inserted directly into the Einstein 
equations. The main result in these works is that black 
holes form with a mass that is bigger than some strictly 
positive value. However, in that approach it is difficult to 
verify whether the implemented corrections break diffeo- 
morphism invariance, which is a fundamental property of 
the theory. 

In contrast to these earlier works, the effective quan- 
tum gravity corrections studied in this letter are via a 
Hamiltonian formulation, and are such that the Dirac 
constraint algebra closes. This means that the diffeo- 
morphism symmetry group is anomaly free. The model 
can be viewed as describing a deformation of the spher- 
ically symmetric sector of the GR plus scalar field the- 
ory, where the deformation comes from a class of possible 
quantum effects. Such an approach may be compared to 
the Ginzburg-Landau theory, where an effective macro- 
scopic description of a quantum theory with some con- 
straints is postulated. 

A numerical study of the resulting field equations re- 
veals qualitative behaviour similar to that found in classi- 
cal GR: there are sub- and supercritical regions of initial 
data characterized by black hole formation and scatter- 
ing, respectively. However, the scaling behaviour in the 
latter region is drastically different, giving a new type of 
oscillation and a mass gap. This is a direct consequence 
of the effective quantum gravity corrections. 

We start from the ADM Hamiltonian formulation of 
GR with the metric given by 

g^^dx^'dx'' = -N^dt^ + A^{N''dt + drf + B'^dVi^, (1) 
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where N is the lapse function and N'^ the radial and sole 
nonzero component of the shift vector. The functions A 
and B define the spatial line element A^dr^ + B'^dD? on 
the topologically fixed space [0, oo) x S"^ . The symmetry 
reduced canonical action of this systems is 

pti POO 

S= / {pAA + pBB+p^<i)~'H) dr dt, (2) 
where the nonzero, equal-time Poisson brackets are 



(3) 



{Ait.,r),pAit,f)} - {Bit,r),pBit,f)} = 
= {<l>{t,r),p^(t,f)} = d{r,f). 

The total Hamiltonian constraint density is 

3 

n = NY,^± +N'^^\r, (4) 

i=l 

where, having chosen the matter Lagrangian density 
^m^_VEKgf^'^0^^0^^^ |g|^det([5]), (5) 

the Hamiltonian and diffeomorphism constraint densities 
respectively are 

n^l^^^iPAA-2pBB) 



- 7^[A'B^' - A{2BB" + + A% 



(6) 
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Here, the part H.'^f' of the Hamiltonian constraint density 
provides intrinsic and extrinsic geometric information of 



[0, oo) X 5^, the part describes the kinetic energy of 
0, and the part its gradient energy. 

We note that 'h!"^' and H^'' are proportional to l/B^, 
because this factor comes from the inverse determinant 
of the spatial metric. The classical fall-off conditions on 
B in the asymptotic region imply that this factor is sin- 
gular at r = 0. However, in the corresponding quantum 
theory the operator has a bounded spectrum 

on the polymer Hilbert space. These considerations to- 
gether with those in ^ motivate the effective constraint 
density 



3 
i=l 



(7) 



where the functions = (5(i)(i?) are designed to ren- 
der the factors 1/B^ regular. The functions can be 
thought of as providing effective quantum gravity cor- 
rections, using the heuristic guidance that operators are 
replaced by expectation values in suitable states. 



Given the modification of the constraints it is essential 
to check their algebra. If the condition 



?(1) - Q{2)Q{3) 



(8) 



is satisfied, the Dirac algebra, that is here the algebra 
generated by the modified constraint densities, closes and 
is anomaly-free. This implies an invariance under diffeo- 
morphisms, that all degrees of freedom are accounted for, 
and that the ADM energy is conserved. 

Since H^^^ is not proportional to 1/B^ the function 
(5(3) can be set equal to 1. To further simplify Hamilton's 
equations of motion generated by H^, the Schwarzschild 
(or polar areal) gauges B = r and pA = can be fixed. 
The line element then takes the familiar form 

g^^dxt'dx" = -N^dt^ + A^dr^ + r^dfl^ (9) 

and the field equations become 



(10) 
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The function 
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encodes the remaining quantum gravity modification and 
s = q^A/{r'^N) is the characteristic speed of the scalar 
field. 

li q — r the modified equations in (jlOp reduce to Ein- 
stein's equations, which are studied in [1] with the result 
that the black hole mass M scales with the amplitude (/)o 
of the initial scalar field profile according to 



M oc \4>o-(t>l\ 



(12) 



In this relation (/)q is such that 0o > gives rise to 
a black hole whereas < (/)g does not. If (/)o = (/)q a 
naked singularity is the result of the collapse. The critical 
exponent 7 ~ 0.374 P, [l3| is universal in the sense that 
it is insensitive to changes of the initial scalar field data. 

For the modified equations, the numerical analysis is 
simplified in coordinates (u, v) adapted to the character- 
istics c^9^ = sdt ± dr of the scalar field. The norm of 
these characteristics is 



3(c±,c±) 



A' 



(13) 



so that lines of constant u, v G R define spacelike di- 
rections if q < r, null directions if q = r, and timelike 
directions if q > r. That is, the scalar field follows null 
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lines of the metric if and only if the modifications are 
switched off. This makes possible a violation of the dom- 
inant energy condition. 

Setting = F/ f the field equations in (|10p can be 
shown to be equivalent to 
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dv4> ^ 



d r = ~- 
2 ' 

= (v^ 

q 



(14) 
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in characteristic coordinates (m, v) (Tlj. If g = r these are 
once again Einstein's equations but this time in double 
null_coordinates. (They have been numerically studied 
in with an economical method.) The difference be- 
tween the characteristic equations in (fT4|) and the effec- 
tive quantum gravity corrected double null equations in 
Q is that here the equations originate from a Hamilto- 
nian formulation. In the remainder of the letter numer- 
ical solutions to the equations in (fT4|) are constructed, 
with the aim of outlining how the classical mass scaling 
law (fT2l) is modified. 
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Figure 1. An illustration of the classical and the modified 
mass scaling law. The relevant features are a mass gap, non- 
sinusoidal oscillations, an overlap with the classical data, and 
a scaling with A. 

At this point the choice of the deformation function q 
remains to be made. It is restricted by the demand that 
Q(i)/r^ be regular at r = and the intuition that quan- 



(1)/ 

tum gravity corrections vanish for r 
choice is therefore q = \/X^ + r 



oo. A possible 
where A is a length 



scale that specifies the extent of the correction. This pa- 
rameter may be thought of as related to the Planck scale 
which is inherent in the polymer quantization method. 
For the following analysis it is helpful to introduce the 
dimensionless variables 



u/A, 



r/A, 



q/X, (15) 



which also render the field equations in (1141) dimension- 
less. The deformation function q is now 



q = vT 



(16) 



The numerical method is based on a uniform grid code 
that uses the unique combination of the Richtmyer two- 
step Lax-Wendroff method and the second-order Runge- 
Kutta method to solve the equations in (fT4|) . That is, 
unlike in previous work, no use of a mesh refinement is 
made. The boundary conditions are r(u, u) = and that 
<i> is regular at v = u. The initial data is r(0, v) = v/2 for 
the geometric sector and 



(/.(0,v) = 00 
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■ exp 



(v-vo)^ 



(17) 



for the matter sector. 

Black hole formation is signalled when the outward 
null expansion 9 = y^2f /F' /r approaches 0. The last 
recorded minimal value of 6 defines a point (u, v) on the 
computational domain, which in turn defines the mass 
M = r(u, v)/2. Figure [1] displays the results of the analy- 
sis. It shows both the classical GR and quantum modified 
mass scaling law (|12|) . The main features of the latter 
are (i) a strictly positive lower bound on the mass M, (ii) 
oscillations, (iii) an overlap with the classical data for 
log]^o((/)o — 4>q) > —4, and (iv) a scaling with the length 
scale A (recall that M = M/X). 

These features are "universal" in the sense that they are 
unaltered if the initial scalar field profile (fT7|) is replaced 
by 



0(0, v) = 00 tanh 



Vo 



(18) 



The value of the mass gap as well as the period and the 
amplitude of the oscillations agrees with the correspond- 
in^ value in Figure [T] up to a relative error of about 1.5% 

Another aspect of the universality is the apparent in- 
dependence of the results on the deformation function q . 
To see this, a one-parameter family of twice differentiable 
functions q can be defined by 



qo + 3(l-2qo)r2 
+ (-3-H8qo)r3 + (l-3qo)r4 



r<l, 
r > 1. 



(19) 



For clarity it is chosen to reduce to the case of GR for 



r > 1 and it is plotted in Figure 2(a) for two sample val- 
ues of qo = q(0) (which must be less than or equal to 1/2 
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Figure 2. Two members of a family of functions q are displayed in Figure 2(a) Fo r r < 1 the parameter qo = 1/4 is at the 
interface between space- and timelike shoes. The mass scaling data in Figure 2(b) displays the effect of a change in qo and 
implies a robustness of the general features. 



for regularity of the equations in (|14p). The mass scal- 



ing relations displayed in Figure 2(b) exhibit the same 
features as in Figure [T] This suggests that the results 
presented here are not only universal with respect to dif- 
ferent initial data profiles but also forminvariant under 
changes of the deformation function. 

In summary, the main developments reported here are 
(i) an application of the effective constraints method 
[2| to spherically symmetric gravitational collapse in the 
ADM formalism motivated by the polymer quantization 
of this sector [6] , and (ii) a numerical simulation of the 
resulting equations. 

The results indicate that the process of black hole for- 
mation is significantly altered if quantum gravity effects 
are taken into account: black holes form with a non-zero 
mass and this mass stays constant up to the observed 
oscillations and until the collapsing matter reaches a cer- 
tain threshold. Beyond this point the classical results are 
recovered. It is interesting that the requirement of reg- 
ularity of the inverse metric function factors, combined 
with closure of the constraint algebra, lead to a univer- 
sality of the results with respect to both the initial data 
and the deformation function. 

It would be interesting to perform a careful analysis of 
the modified critical solution. Classically, this solution is 
a naked singularity. In the effective theory it is possible 
that this is replaced by an unstable "boson star" solution 
which sits at the threshold of the supercritical regime. 
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